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I In the present paper we outline the stochastic Umit approach to superfluidity. The 

^ ■ Hamiltonian describing the interaction between the Bose condensate and the normal phase 

Tlj" . is introduced. Sufficient in the stochastic limit condition of superfluidity is proposed. 

I Existence of superfluidity in the stochastic limit of this system is proved and the non- 

' linear (quadratic) equation of motion describing the superfluid liquid is obtained. 

o : 
^ ■ 

^ ■ 1 Introduction 

^1 The theory of superfluidity was developed by L.D. Landau and N.N.Bogoliubov, cf. [|1|, 0, 
I 0) HI) for an introduction see p. An analogue of the approach of 0, 0] was applied to 
^ ■ superconductivity theory, cf. p. 

O^. The essence of the superfluidity phenomenon is that the Bose condensate becomes superfluid 

J> ! and friction between condensate and normal phase disappears. N.N.Bogoliubov in [^j, |Q found 
I that this can be explained as an effect of stabilization of the condensate by interaction between 
^ [ particles. 

In the present paper, using as a starting point the approach of 0, 0, we introduce an 
Hamiltonian that describes the interaction between the Bose condensate and the normal phase 
and investigate this Hamiltonian using the stochastic limit approach cf. |^. 

The name stochastic limit is due to the property that in this approach the quantum field 
is approximated by a quantum white noise and the Schrodinger equation is approximated by a 
white noise Hamiltonian equation. In the stochastic limit we start from an Hamiltonian of the 
form 

H = Ho + XHi 

where A is a coupling constant and we make the time rescaling t ^ t/X^ in the solution of the 
Schrodinger equation in interaction picture U^'^^ = e**^oe~**^, associated to the Hamiltonian 
H, i.e. 
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with Hj(t) = e^^^°Hie (the evolved interaction Hamiltonian) . This leads to the rescaled 
equation 

and one wants to study the hmits, in a topology to be specified, 

\imHj{t/X'^) = Ht 

Moreover one wants to prove that Ut is the solution of the white noise Hamiltonian equation 

dtUt = -iHtUt ; t/o = 1 

which is equvalent to a quantum stochastic differential equation. 

The structure of the present paper is as follows. In section 2 we discuss the Bogoliubov- 
Landau condition of superfiuidity. In sections 3 and 4 we review some basic facts from |^ on 
Bose condensation in ideal and non-ideal Rose gases. In section 5 we propose the Hamiltonian 
of condensate-normal phase interaction and investigate this model using the stochastic limit. 
In section 6 we construct by the stochastic golden rule the master equation for the condensate- 
normal phase interaction and prove the existence of superfiuidity in this system. We prove that 
the standard Bogoliubov-Landau condition is not sufficient to garantee the superfiuidity in the 
stochastic limit and introduce a procedure to overcome this problem. In section 7 we sum up 
the main conclusion that can be draw from our results. 



2 The Landau and Bogoliubov ideas on superfluidity 
theory 

Tissa and London conjectured that the existence of the Bose condensate can explain super- 
fiuidity in the sense that the condensate corresponds to the superfiuid component and the 
temperature phase corresponds to normal component. 

By Landau's argument the existence of a condensate itself is not enough to prove 
superfiuidity and moreover there is no superfiuidity in ideal Bose gas. 

Consider the condensate moving with velocity u and the non-moving normal state. Suppose 
that the friction of the condensate with the normal phase stimulates the transition of a particle 
in condensate from the state with velocity u to the state with velocity u — —. Then this 

particle will have energy ^"'^"^^ . Before the transition the particle had energy Since the 
transition is due to friction with normal state this corresponds to the excitation of the particle 

2 

of the normal phase with energy We get the energy difference 

Energy after transition + Energy of excitation - Energy before transition = 

imu — pY mu^ 

= 1 = pu < for small p 

2m 2m 2 m 
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This means that such transitions are energetically possible and that the condensate is non- 
stable. 

Bogoliubov found that the interaction between the particles in thecondensate may 

stabilize the condensate and modify the energy difference in such a way that transitions to 
slower motion due to friction will proceed with positive energy difference and therefore will be 
energetically forbidden. 

Bogoliubov found that interaction makes the state of normal phase instable. To stabilize 
the state one has to make a Bogoliubov transformation. This procedure changes the dispersion 
E{k) of excitations of normal phase, cf. [Q, 0], 0, 0. 

After the canonical transformation the energy difference for transition p p — k is given 

by 

E{k) + e{p - k) - sip) (1) 

2 

where e{p) = ^ and E(p) has different dependence on p. In this case one can overcome the 
Landau objection. Consider the following examples. 

EXAMPLE 1: Consider radiative dispersion E{k) = c\k\ + . . .. In this case E{k) for small 
enough k is proportional to \k\ and takes the form 

c\k\ - pk + > 0, \p\<c (2) 
EXAMPLE 2: Consider the polaron model with E{k) = uj. In this case (|l|) becomes 

uo -kp+ > 0, \p\ < y/uj 

Let us show that the results of Landau and Bogoliubov discussed above are connected with 
the stochastic limit approach. One of important properties of the stochastic limit is that the 
main properties of the dynamics is controlled by 5-functions of energy differences of the form 

6{E{k) + e{p - k) - e{p)) (3) 

(this coincides with the 6 of (|I|)). The possibility to control the dynamics by exploiting such 
5-functions was called in ||^ the Cheshire Cat effect. We find that the Bogoliubov-Landau 
condition of superfluidity is connected with the Cheshire Cat effect developed in the stochastic 
limit approach. 



3 Condensation of an ideal Bose gas 

In the present section we discuss standard material on Bose condensation of an ideal Bose gas, 
cf. [^. The Hamiltonian of an ideal Bose gas in second quantization is given by 

H = J uj{p)b*{p)b{p)dp 

The grand canonical ensemble Hamiltonian is 

T = H-i^N, N = I b*{p)b{p)dp 
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The equilibrium state of the gas is the unique mean zero gauge invariant state with expectations 

{b*{p)bip)) = {n{p)) = ^,(,(^)^^, _ ^ 

giving the density of particles with momentum p and where /5 = | is the inverse temperature. 
We take /i < since {n{p)) > 0. 

The integral over the density of degrees of freedom of gas is equal to the density of the gas: 

{n{p))dp = p 

To describe Bose condensation, we consider the zero mode of the field 69 as singled out from 
the other modes 

h,K] = i 

and take, as the reference state of the field, the state not the equilibrium one but a gaussian 
state of the form ^ 

{nip)) = c5{p) + _ ^ 

Notice that the chemical potential in this formula is taken equal to zero. 

We consider the constant c = c{9) depending on the temperature. To determine this 
dependence we use the conservation of number of particles, given by the following: 

^ dp = p (4) 



3/3w(p) _ 1 

Therefore 

/ \ 2 



„2 



where c? = 3 and uj{p) = 6 < Oc- The critical temperature 9c is the temperature when the 
integral in @) becomes equal to the total density of particles p and condensate disappears. 
This implies the state of the Bose gas in the form 

{<v)) = p\l-\^\\5{p)+ ^ 



4 Condensation of a non— ideal Bose gas 

In the present section we review material from containing the discussion of the Bose conden- 
sation of non-ideal Bose gas and prove the applicability of the Bogoliubov-Landau condition 
of stability of the condensate. 

The second quantized Hamiltonian of a non-ideal Bose gas (for the grand canonical ensem- 
ble) is: 

r = - fJ-)b*ip)b{p) + ^ J2 9i.Pi - P'i)b*ipi)b*{p2)b{p2)b{p[) 

P Pl+P2=p'i+P2 
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We consider a state such that almost all particles with p = are in the condensate. This imphes 

(6*60) = No, No » 1 
where we denote the creation of the particle with p — 0. Since 

(bob*,) = No + l 

and A'o + 1 is almost equal to Nq, we can consider 60 as a classical (commuting) variable with 

l&or = blbo = b*obo = No, bo = ^/N'oe"^ 

This allows to simplify the Hamiltonian P. Since b{p) with p ^ are not in the condensate, 
amplitudes of b{p) are small with respect to bo- Therefore we can keep in the Hamiltonan F 
only the terms of second order in b{p), p ^ and skip higher order terms (3 and 4). 

This mean field approximation will be valid only when almost all particles are in the con- 
densate, i.e. for a temperature close to zero: 9 = 0. 

In this approximation the fourth order term in F takes the form 



A f ^(o)7V„2 + Y^g^p)b(p)b{-p) + blbo + 9i-pmp)bip)+ 

+2b*oboj:9(0)b(p)b(p) + blJ2g(p)b*(p)b*(-p) 



Using the conditions g{p) = g{—p), a;(0) = we get for F 

r = ^9{0) - uNo + E (u{p) - + ^^(0) + ^9{p?j b*{p)b{p) + 

+ ^K'E9iP)Kp)K-p) + ^blY,g{p)b*{p)b%-p) 
This imphes that the ground state of F has the energy 

El = ^^(0) - f^No 
The chemical potential is defined by the equation 

dNo 

and in our case it is equal to 

^^^-^9{'^) 

With this chemical potential and making the canonical transformation b{p) b{p)e^^ we cancel 
the phase of the complex number 60 — \/Noe.^'^ so that the Hamiltonian F becomes 



r = -^5(0) + E + ^9{P)) h*{p)b{p) + ^ ^M{b*{py{-P) + K-p)Kp)) 
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Let us make a canonical transformation to diagonalize this quadratic Hamiltonian. We consider 
the following transformation 

a{p) = Upb{p) + Vpb*{-p) 
a*{p) = Upb*{p) + Vpb{-p) 

- = 1) = W-P) = ^-p; Wp) Vp eK 
The inverse transformation is given by 

Kp) = Upa{p) - Vpa*{-p) 

b*{p) = Upa*{p) - Vpa{-p) 

After this canonical transformation the off-diagonal terms in the Hamiltonian (coefficients 
of a*{p)a*{—p) and a{—p)a{p)) are equal to 

- UpVp (^{p) + ^^(P)^ + + vl)^^g(j)) (5) 

and the diagonal terms (coefficients of a*{j))a{p)) are equal to 

[Up + Vp) \ uj{p) + -^g[P) 1 - ^UpVp^gip) (6) 

Since Up — = 1 we can use the hyperbolic parametrization 

Up = ch X, Vp = sh X 

The compensation equation (vanishing of off-diagonal terms (§)) in this parametrization takes 
the form 



th 2x = ——^ 



iuip) + ^g{p) 



This implies 



Up + Vp= ch 2x 



^^0 I \ 



^9{p) 

J-pt'p — Dli ^U. — . _ 



2UnVr, = sh 2x 



V 



'1- 



Therefore for (BI) we get 



Eip) = \Jco\p) + ^u;{p)g{p) (7) 

Condition g{0) > that provides positivity of the value under the square root corresponds to 
the domination of repulsion, cf. 0. 
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The Hamiltonian takes the form 



r = -^9iO) + Y.E{pX{p)a{p) (8) 

The state for hte new Hamiltonian @ corresponds to a state of the non-ideal Bose gas 
modified due to interaction. 

Let us discuss the Bogoliubov-Landau condition for Hamiltonian the (|^). We get for (|1|) 

m (u — —) mip rp- 
Eip) + ^-^^-^ = E{p)-pu + ^ (9) 

Since E{p) = c\p\ for small p we can get the situation when for u < c the energy of excitation is 
positive and transition from condensate to normal state proceeds with consumption of energy. 
In this case we can apply the discussion of (0) which shows that the Bogoliubov-Landau 
condition is applicable to (H). 



5 Condensate— normal phase interaction in the stochas- 
tic hmit 

In the present section we will discuss how to describe the Bogoliubov-Landau condition by the 
stochastic limit approach. 

We will consider the stochastic limit of a system with phase transition. We will see that the 
phase transition leads to arising of non-linear master and kinetic equations. This shows that 
phase transition and the stochastic limit procedure do not commute: to investigate a system 
with phase transition by the stochastic limit we should describe first the phase transition. 

We introduce the Hamiltonian for condensate-normal state interaction and investigate it in 
the stochastic limit. 

Consider a system with two Bose fields (phases) 

[c(p), c+ip')] = 5{p - p'), [aik), a+ik')] = 5{k - k') 

We describe the condensate by the Bose field c(p) (system degrees of freedom) labeled by a 
velocity index p with free Hamiltonian 

Hc = j e{p)c*{p)c{p)d^p 

with dispersion 




and state 

{c*ipW)) = N{p)5{p-p') 
This state describes the distribution of condensate over velocities. 



7 



The second Bose field a{p) (reservoir degrees of freedom) describes excitations of the non- 
ideal Bose gas (normal phase) which diagonalize the interaction (considered in the previous 
section). This excitations correspond to pairs of particles with opposite momenta and has free 
Hamiltonian 

Hns = J E{p)a*{p)a{p)d^p 
The state of a{p) is an equilibrium 

The total Hamiltonian will be 



H = H, + H^s + XHi (10) 

with the interaction 



Hi = J J f{k,p)c*{p)c{p~ k)a{k)d'pd'k+ h.c. (11) 

where f{u,p) is the form-factor (complex valued test function). 

After this we apply to the Hamiltonian ( [T0| ) the stochastic limit approach: by the stochastic 
golden rule we construct the master equation for the density of particles in the condensate 
c*{p)c{p) and prove that the stochastic golden rule gives the Bogoliubov-Landau condition for 
this system. 

This approach is close to the discussion in Abrikosov-Gorkov-Dzialoshinskii book: Consider 
the Bose liquid with zero temperature that flows with velocity v. In the presence of friction in 
the liquid there will be elementary excitations with different velocities. An analysis of these 
elementary excitations gives the Bogoliubov-Landau condition, cf 0, P, 0. 

The friction is due to the interaction term (|ll]). The rescaled free evolution of the interaction 
(|1T]) (of the collective field) is equal to 

Ax{p, Kt) = \ e-^c^{p)a{k)c{p - k)e-'^ = \ c-^{p)a{k)c{p - A;)e-^*^(P'^-)/^' (12) 
A A 

A^ip, k,t) = j e'^c+ip - k)a+{k)c{p)e'-^ = i c+(p - A;)a+(A;)c(p)e'*^(^''=)/^' (13) 
where 

E{p,k)=E{k)+e{p-k)-e{p) (14) 

is the corresponding energy difference. 

The difference of energies (|l^) coincides with the difference of energies in the Bogoliubov- 
Landau condition (|^) that justify the fact that Hamiltonian ( [T0| ) describe the non-ideal Bose 
gas considered above. 

The stochastic limit of the rescaled evolution of the interaction (|I^) gives rise to a quantum 
white noise (master field): 

lim A\{p, k, t) = B{p, k, t) 

A — 5-0 
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After the stochastic hmit we cannot separate in the master field the degrees of freedom of 
reservoir and normal state. This means that the degrees of freedom of the condensate and of 



the normal phase become entangled even at kinematical level, cf. . 

The master field B{p, k, t) will satisfy the variant of quantum Boltzmann statistics given by 
the following theorem 

Theorem 1. The entangled noise algebra is generated by Bi{p, k, t), B2{p, k, t), n{p) where 

B{p,k,t) = B,{p,k,t) + B;{p,k,t) 

with relations 

Bib; = B2BI = 

B,{p, k, t)Bl{p', k', t') = 27cSit - t')S{E{p, k))n{p){N{p - k) + - p') ^^^~^X) ^^^^ 

B2{p,k,t)B;{p',k',t') = 27T6{t-t')6{E{p,k))n{p-k){N{p) + l)6{p-p')^^^ (16) 

[n{p'), Bfip, k, t)] = {±){6{p' -p)- 6{p' -p + k))Bf{p, k, t) 
[n(p'), B^{p, k, t)] = iTmp' -p)- 5{p' -p + k))B^ip, k, t) 

[n{p),n{p')] = 

Here Br = Bi, 5+ = B*. 



The zero temperature version of this theorem was obtained in |10|. The one particle sector 



relations {n{p) = 1) were investigated in |jTT|, |T^, [|T3| 



Consider now the evolution in the stochastic limit for the considered Hamiltonian. By the 
stochastic golden rule, cf. [0] we get 

Theorem 2. The stochastic differential equation for Ut has the form 

dUt = (^-i J dkdp(j{k,p)dB{p,k,t) + dB\p,k,t)f{k,p)) - 

-dt{f\f)_-dtJJ\r)+)Ut (17) 
where the stochastic diferentials satisfy the Ito table 

dB{p, k, t)dB\p\ k', t) = 27m{p)S{uj{k) + e{p - k) - e{p)) 

{N{p -k) + mp - p')-^^^^^dt (18) 
dB^{p, k, t)dB{p\ k', t) = 27in{p - k)5{uj{k) + e{p - k) - e{p)) 

{N{p) + mp-p')^^^dt (19) 



(/!/)_= dkdp\f{k,p)\ 



2 -in{p){N{p- k) + 1) 



uj{k) + e{p - A;) - e{p) - iO 1 - e-Z^^C^) 
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uj{k) + e{p -k) - e{p) - iO e^^C^) 



One particle version of this quantum stochastic differential equation was obtained and in- 
vestigated in 



6 Master equation and superfluidity 

In the present section we consider the master equation for the stochastic dynamics described in 
the previous section. The master equation is the equation for the expectation over the degrees 
of freedom of the reservoir (normal state in the considered case) of the Heisenberg dynamics 
for some observable X 

1^' - " - - (20) 

with Ut defined by theorem 2. 

Consider the master equation for the density n{p) = c*{p)c{p). The stochastic golden rule 
(application of theorems 2 and 1) gives 

j^ritiq) = J dkdp\f{k,p)\\6{q - p) - 6{q - p + k))2rt6{E{k) + e{p - k) - e{p)) 

(n,{p - k){N{p) + - n,{p){N{p - k) + 1)^^^^) (21) 

Since nt{p) is in the abelian subalgebra this equation is classical and we can consider n{p) as a 
classical distribution. 

Taking the integral over p one gets 

d 



-ntiq) = 2n J dk{\f{k,q\^5{E{k)+e{q~k)-e{q)) 
nt{q - k){Niq) + I) _ ^ " ^.(g)(iV(g - k) + 1)^3^^) ^ 



-\f{k,q + k)\'6{E{k)+e{q)-e{q + k)) 
ntiq)iN{q + k) + - nM + k){N{q) + 1)^^^^ ) ) (22) 



Let us discuss the connection of the Bogoliubov-Landau condition with (^). Since the disper- 

2m 



sion (1^ for small \k\ is radiative E{k) = c\k\ and e{j)) = ^ "we get 



I k"^ 1 

E{k) + eiq - k) - eiq) = c\k\ qk H > 0, — Igl < c 

m 2m m 

1 k"^ 

E(k) + e(q) — eiq + k) = c\k\ qk 

m 2m 
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This implies that the first (5-function in ( p2D vanishes for ^ | g | < c. Let us note that the second 
5-function in (^2]) does not vanish for any q. 

To get a superfluid motion in (|2^) it is sufficient to consider such n{q) that the following 
products will vanish for arbitrary q and k: 

n{q)S{E{k) + eiq ~ k) ~ e{q)) = (23) 

n{q - k)6{E{k) + e{q - k) - e{q)) = (24) 
If n{q) would satisfy (p3|), (p4[) then the RHS of (E^) would vanish. 



To satisfy (^) it is sufficient to take the support of n{q) is concentrated in the ball \q\ < mc. 
But (^4]) can not be satisfied when the support of n{q) is concentrated in the ball |g| < mc. 
This implies that (|2^) does not describe a superfluid liquid and to investigate superfluidity in 
the stochastic limit one needs some additional arguments. 

To overcome this problem we propose the following construction. 

First, since N{p) >> 1 we can identify and + 1. 

Second (and most important), we identify nt{p) with expectation {n{p)) = N{p) in (^) (we 
substitute N{p) by nt{p) in this formula). The identification 

{n{p)) = Nip) := ntip) (25) 

means that we consider the stochastic limit with the state depending on time in the slow time 
scale of the stochastic limit. This modifies the stochastic limit procedure in order to take into 
account the effects of phase transition. 

This condition may be justified by the following argument. Equation ( pO|) for X = n{q) is 

d_ 

dt 



(UMqn) = {U9o{n{qm (26) 



where 9o{n{q)) is equal to the RHS of (|^). Then we apply the identification (|25| ) for t = to 
(pBI). We get, in the integral on RHS of (p6D, combinations of terms of the folowing form 



it{n{q)n{q - k)) = ]t{n{q))]t{n{q - k)) 
since jt is a homomorphism. After this we apply the semiclassical approximation 

{jt{n{q))jt{n{q - k))) = {jt{n{q))) {jt{n{q - k))) = nt{q)nt{q - k) 



that justifies ( |25D for any t. 

Condition (pSj) introduces the self-interaction into the considered model and makes equation 
( p2D nonlinear. Condition (|25| ) implies for ( [2^ ) the following: 



d 
di 



nt{q) = -2n J dk{\f{k, q\H{E{k) + e{q - k) - e{q))nt{q)nt{q - k)- 

-\f{k, q + k)\'6{E{k) + e{q) - e{q + k))nt{q)nt{q + A;)) (27) 



Let us note that this equation is nonlinear (quadratic). 
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For ( p7D the condition of superfluidity is reduced to (^) (condition (^) can be ignored). 
This imphes that (pT]) describes a superfluid flow. We get that in the stochastic limit the 
condition 

Supp n{q) C {q : \q\ < mc] (28) 

is sufficient for superfluidity in the sence that, under this condition, the right hand side of ( pT]) 
is zero and the density nt{q), of the condensate, is constant. 

We call equation ( pTj) the equation of motion of superfluid liquid. 

Non-linear master equation ( ^7]) is an example of a general phenomenon. Application of the 
stochastic limit to systems with phase transitions will generally create non-linear master and 
kinetic equations. Non-linearity will enter through the self-interaction given by an analogue 
of ( P5| ) for a certain phase (an analogue of the condensate considered in the present paper). 

7 Conclusion 

We summarize our main results as follows: 

1) The Bogoliubov-Landau condition of superfluidity follows from the stochastic limit approach. 

2) The stochastic limit approach gives not only a condition of superfluidity but also a natural 
candidate for the equation of motion of the superfluid liquid (pT)). 

3) This equation is nonlinear (quadratic). 

4) Without the introduction of a non-linearity (identification of n{p) and N{p)) it is impossible 
to get superfiuidity in the stochastic limit approach. 

5) Nonlinear master equations should be a general feature of the stochastic limit of self- 
interacting systems. 
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